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A  diapycnal  mixing  budget  on  the  Oregon  shelf 

Eric  A.  D’Asaro^ 

Applied  Physics  Laboratory  and  Department  of  Oceanography,  1013  NE  40th  Street,  Seattle,  Washington  98105 
Abstract 

Although  isopycnal  mixing  is  undoubtedly  important  at  global  and  gyre  scales,  the  relative  importance  of 
isopycnal  and  diapycnal  mixing  on  much  smaller  scales  is  uncertain.  This  issue  is  investigated  using  35  d  of  data 
from  a  Lagrangian  float  deployed  on  a  mid-depth  isopycnal  on  the  Oregon  shelf.  Measurements  of  temperature, 
salinity,  and  pressure  maintain  the  float  on  the  isopycnal;  its  high-frequency  diapycnal  deviations  are  used  to 
estimate  the  diapycnal  diffusivity  using  an  inertial  subrange  method;  lower-frequency  deviations,  including 
intentional  profiles  to  the  surface,  are  used  to  estimate  diapycnal  derivatives  near  the  target  isopycnal.  Downward 
irradiance  at  490  nm  is  used  to  calibrate  chlorophyll  fluorescence  measurements  and  compute  solar  heating  rates. 
Corrections  for  the  diapycnal  deviations  provide  a  nearly  continuous  isopycnal  time  series  of  spice  (a  temperature 
and  salinity  combination  nearly  orthogonal  to  potential  density)  and  chlorophyll.  A  new  formulation  of  the 
diffusion  equation  in  isopycnal  coordinates  is  derived  and  used  to  test  the  accuracy  of  purely  diapycnal  mixing 
balances  for  spice  and  chlorophyll.  On  vertical  scales  of  about  10  m  and  timescales  of  about  2  d,  isopycnal  spice 
variations  are  mostly  controlled  by  diapycnal  mixing,  although  other  processes,  presumably  isopycnal  mixing,  are 
sometimes  important.  Processes  other  than  diapycnal  mixing  control  isopycnal  chlorophyll  variations  on  these 
scales.  Likely  candidates  include  isopycnal  mixing  with  a  nearby  bloom,  planktonic  sinking  out  of  this  bloom,  or 
possibly  local  phytoplankton  growth.  Thus  both  isopycnal  and  diapycnal  mixing  can  be  important  at  these  small 
scales. 


Isopycnal  mixing  is  clearly  important  to  the  large  and 
mesoscale  distribution  of  scalars  in  the  ocean.  Gyre 
circulations  and  mesoscale  eddies  stir  the  global-scale 
gradients  in  temperature  and  salinity  to  smaller  scales.  On 
these  large  scales  the  flow  is  nearly  adiabatic  and  the 
stirring  is  therefore  nearly  isopycnal.  Both  models  and  data 
suggest  that  such  stirring  moves  scalar  variance  down  to  at 
least  kilometer  scales  (Sundermeyer  and  Price  1998). 
Diabatic  mixing  removes  this  variance  on  smaller  scales. 
In  models,  this  is  often  done  by  isopycnal  or  horizontal 
diffusivity,  but  it  is  not  known  whether  isopycnal  or 
diapycnal  mixing  dominates  in  the  real  ocean.  Process 
studies  of  ocean  mixing  have  concentrated  almost  exclu¬ 
sively  on  diapycnal  mixing  at  scales  of  meters  to 
millimeters.  Established  techniques  for  the  estimation  of 
diapycnal  mixing  using  microstructure  measurements  rely 
on  diapycnal  (or  vertical)  variance  balances  (Gregg  1987; 
Thorpe  2005)  with  little  reference  to  the  possibility  of 
isopycnal  processes  playing  an  important  role.  One 
possibility  is  that  the  flux  of  scalar  variance  from  large  to 
small  scale  occurs  isopycnally  at  large  scales  and  diapyc- 
nally  at  microscales,  with  the  transition  scale  being  in  the 
poorly  measured  range  between  meters  and  kilometers.  If 
so,  then  the  local  Lagrangian  rates  of  change,  i.e.,  the 
turbulent  flux  divergences,  would  be  balanced  diapycnally 
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when  viewed  at  meter  scales,  but  isopycnally  when 
averaged  over  larger  scales.  In  this  paper,  the  relative 
importance  of  isopycnal  and  diapycnal  mixing  at  meter 
scales  is  tested  by  comparing  the  Lagrangian  rate  of  change 
of  the  scalar  “spice,”  a  combination  of  temperature  and 
salinity  that  represents  water  mass  properties,  with  the 
divergence  of  the  diapycnal  flux  of  spice.  A  similar  test  is 
made  for  chlorophyll  concentration.  The  following  sections 
describe  the  theoretical  approach,  the  instrumentation  and 
data  processing,  the  tests,  and  their  interpretation. 

Techniques 

Mixing  rates  from  an  isopycnal  budget — McDougalTs 
analysis:  Consider  an  “isopycnal”  float,  i.e.,  a  float  that 
moves  horizontally  with  the  velocity  of  the  water  and 
moves  vertically  to  remain  on  a  surface  of  constant 
potential  density  referenced  to  the  surface,  g{S,T,P),  with 
value  (To-  The  equation  for  the  rate  of  change  of  a  scalar  C 
measured  on  the  float  is 

DC 

=  C,  =  -{e  -  Vd)Cd  +  VCm  +  Vi-KViC  +  Q%\) 
D; 

where  the  subscript  d  implies  a  diapycnal  partial  derivative, 
V/  is  an  isopycnal  horizontal  graciient  operator,  V  is  the 
diapycnal  diffusivity,  K  is  the  isopycnal  diffusivity,  and  e  is 
the  diapycnal  velocity,  i.e.,  the  velocity  across  or  through 
the  isopycnal  surface.  Other  sources  of  C,  for  example, 
solar  heating  for  a  temperature  budget,  are  represented  by 
The  actions  of  e  and  Vd  cannot  be  distinguished  in 
budget  calculations.  Accordingly,  an  effective  diapycnal 
advective  velocity  is  defined  as 

e  =  e  -  Vd  (2) 
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McDougall  (1984)  takes  these  equations  for  6  and  S  with 
=  0  and  =  0,  eliminates  e,  and  uses  the  relationship 
(his  equation  6) 


where  Q,  with  no  superscript,  is  the  local  rate  of  heating  in 
units  of  6t  [°C  s“i].  Because  the  measurements  are  along  an 
isopycnal, 


Soo  = 

9ddSd(9dSdd  \ 

9d  \^d  9dd  ) 

(3) 

II 

o 

(8) 

to  derive 

so  that 

0,  =  VgN~ 

-^9]pSee  +  K{V]9  +  N.L.) 

(4) 

d^^^dd  Q 

e  =  U —  —  aa  — 

<^d  O'rf 

(9) 

relating  the  rate  of  change  of  temperature  measured  by  an 
isopycnal  float  to  the  diffusivities  and  thermohaline 
structure.  Here  P  is  the  haline  contraction  coefficient  P  = 
P~^8sp\o.p  and  N  is  the  buoyancy  frequency. 

In  general,  both  the  diapycnal  and  isopycnal  diffusivity 
can  change  the  temperature.  Isopycnal  mixing  also  acts 
through  the  nonlinearity  of  the  equation  of  state  to 
contribute  the  terms  labeled  “N.L.”  These  are  described 
more  fully  in  McDougall  (1984).  For  K  =  0,  there  is  a  direct 
relationship  between  6,  and  V,  which  depends  on  the 
curvature  of  the  T/S  relationship  Soo-  Practically,  however, 
there  are  problems  with  Eq.  4.  Its  form  is  not  intuitive  and 
Soo  can  be  zero,  infinite,  or  even  multivalued.  Accordingly, 
an  alternative  form  is  sought. 

The  spice  variable  H:  Potential  density  can  be  expanded 
locally  in  terms  of  p  and  a,  the  thermal  expansion 
coefficient,  as 

a{S,0)  =  u(So,0o)  -  a(d  -  do)  +  P{S  -  So)  (5) 

Flament  (2002)  defines  spiciness  (abbreviated  here  as 
spice),  denoted  by  n(S,  9),  as  the  natural  companion  to 
potential  density  (j{S,  9).  The  units  of  u  and  H  are  the  same 
and  n  is  defined  to  be  zero  at  0  =  0,  salinity  (S)  =  35. 
Locally,  H  is  expanded  as 


requiring  curvature  of  the  density  profile  or  local  heating  or 
cooling  to  produce  a  diapycnal  velocity.  The  equivalent 
equation  for  H  is 

=  —elld  +  DlVdd  +  WoQ  (10) 

A  diagnostic  equation:  Equations  9  and  10  can  be  solved 
for  V> 

n,  =  VUmU  -  +  alne  +  (11) 

V  Od  \  (^dj 

The  right  two  terms  are  the  usual  diffusion  equation  with  a 
source,  but  with  the  addition  of  the  extra  factors  in 
parentheses  resulting  from  diapycnal  advection  due  to 
diapycnal  diffusion.  Diffusion  moves  the  isopycnal  through 
the  water.  These  extra  factors  quantify  this  effect. 

For  this  analysis,  it  will  be  easier  to  take  a  derivatives 
than  diapycnal  derivatives.  Therefore  H  =  n((T)  will  be 
assumed  for  the  purpose  of  evaluating  the  right  side  of 
Eq.  11.  Obviously  this  cannot  be  exactly  true  if  H,  is 
nonzero  on  an  isopycnal,  so  it  cannot  be  assumed 
otherwise.  Using  the  chain  rule  to  express  Hj  and  I\dd  in 
terms  of  a  derivatives  yields 


n(5,0)  =  n(5o,0o)  +  no(5o,5o) 


(0 


0o)  +-{S 
0( 


(6) 


where  Hq  is  the  partial  derivative  of  H  with  respect  to  9  and 
the  defining  relationship  pHo  =  aH^’  has  been  used.  The 
exact  definition  within  these  constraints  is  formed  requiring 
n  to  be  minimally  correlated  with  a  over  the  range  of 
oceanic  water  masses. 

In  general,  H  is  not  a  linear  function  of  6  and  S  because  a 
and  P  are  functions  of  9  and  S.  Therefore,  like  oo.  U  is,  in 
general,  not  conserved  during  mixing.  These  effects  will  be 
ignored  in  this  paper  and  the  full  nonlinear,  nonconserva¬ 
tive  form  of  the  Flament  (2002)  H  will  be  used.  However, 
parallel  computations  using  a  linearized,  conservative  form 
of  H  were  also  made,  as  described  in  the  section  “Mixing 
budgets.”  The  differences  are  negligible. 


An  isopycnal  equation  for  H:  If  depth  and  temperature 
deviations  are  small,  the  deviations  of  the  equation  of  state 
from  Eqs.  5  and  6  should  be  small.  Applying  Eq.  1  to  0  and 
S  with  K  =  0  and  using  Eq.  5  yields  an  equation  for  a 

^  =  -e(Jd  +  Vffdd  -  acrQ  (7) 


Hr  =  +  Q^Ylo  +  aG\l„) 


or  in  more  conventional  notation 


DU 

~Dt 


gl 


dH 


(12) 


(13) 


This  can  be  used  to  estimate  the  diapycnal  diffusivity 


On 


tlr 


(14) 


where  Vn  is  the  diffusivity  computed  using  Eq.  14  and 


n,  =  H,  -  Q^Yle  +  auH^)  (15) 

is  H,  corrected  for  solar  heating.  These  expressions  contain 
the  same  physics  as  Eq.  4  with  the  addition  of  the  heating 
term,  but  are  written  in  terms  of  diapycnal  derivatives 
instead  of  0  and  S  derivatives.  The  change  in  properties 
along  the  isopycnal  due  to  diapycnal  mixing  depends  on  the 
curvature  of  the  H-u  diagram  and  the  stratification. 
However,  the  form  is  much  simpler,  and  thus  easier  to 
evaluate,  and  the  curvature  is  always  well  defined. 


Diapycnal  mixing  budget 
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Scalar  changes;  A  similar  analysis  applies  to  the  rate  of 
change  of  any  scalar.  From  Eq.  1  using  Eq.  9: 


spectra  in  D’Asaro  and  Lien  (2007)  collapse  to  a  single 
nondimensional  form 


C,  =  V 


^  ^dd  ^ 

Wrf  —  — w 


+  2' 


(16)  =  P,xH{m/o)L)  =  P/X. 


Again  assuming  C  =  C((t)  and  using  the  chain  rule 
Ct  =  VCaa(y]i  +  ao-QQ  +  Q 


(17) 


which,  not  surprisingly,  is  a  slightly  more  general  form  of 
Eq.  12. 


1  + 


3.0(U/^ 


4.5caz, 


(20) 


where  H  is  an  empirical  nondimensional  function,  increas¬ 
ing  at  high  frequencies  because  of  advection  of  density 
fluctuations  past  the  float’s  density  sensor.  Here 


col  =  (e/E^) 


2\1/3 


(21) 


Budget  method  summary:  Equation  12  or  17  describes 
the  local  budget  for  scalars  fl  or  C  in  an  isopycnal 
coordinate  system  moving  with  the  water.  Because  water 
moves  across  isopycnals  in  the  presence  of  mixing,  this  is 
different  from  a  fully  Lagrangian  coordinate  system 
because  of  the  diapycnal  velocity  e.  The  value  of  V  can 
be  estimated  by  these  equations  assuming  no  isopycnal 
mixing.  The  diapycnal  velocity  e  is  implicit  in  the 
calculation,  but  can  be  evaluated  from  Eq.  9  if  the  gradient 
of  stratiflcation  is  known.  Here,  the  strategy  is  to 

compare  the  value  of  t)  estimated  in  this  way  with 
microscale  estimates  and  thereby  test  the  assumption  that 
K=Q. 


Diapycnal  mixing  from  microscale  measurements — The 
Osborn  and  Cox  (1972)  method  computes  the  down- 
gradient  flux  of  any  scalar  from  its  dissipation  rate  and 
its  gradient.  Winters  and  D’Asaro  (1996)  show  that  this 
formula  is  exact  if  the  correct  definition  of  diapycnal 
gradient  is  used.  Evaluation  of  the  difference  between  their 
definition  and  the  usual  spatial  gradient  using  density  fields 
from  numerical  simulations  (D’Asaro  et  al.  2004)  and 
ocean  data  shows  that  there  is  little  difference.  Applying 
this  to  potential  density  a,  so  as  to  yield  the  diapycnal 
diffusivity  of  potential  density,  yields 


2(1)2  ^2AV2 


(18) 


The  dissipation  rate  of  potential  density  variance  x  was 
estimated  as  in  D’Asaro  and  Lien  (2007).  In  a  turbulent 
fluid  Lagrangian  frequency  spectra  of  scalars  d>„(a))  exhibit 
an  inertial  subrange  for  Lagrangian  frequencies  co  above 
about  N.  This  corresponds  to  the  -  5/3  spectrum  found  for 
scalar  wave  number  spectra  (Sreenivasan  1996),  but  here 
found  for  frequency  instead  of  wave  number.  D’Asaro  and 
Lien  (2007)  show  that  this  spectrum  should  have  the  form 


is  the  frequency  at  which  the  float  size  L  becomes 
important.  The  value  of  x  is  estimated  from  the  average 
value  of  0D,{cci) / P^H{co / col)  for  m  <  loj^  (D’Asaro  and 
Lien  2007). 

The  dissipation  rate  of  kinetic  energy  e,  needed  in  Eq.  21, 
was  evaluated  by  equating  Eq.  18  with  the  Osborn  (1980) 
formula 


U  =  r 


e 

w 


(22) 


relating  the  diapycnal  diffusivity  to  £  in  terms  of  a  “mixing 
efficiency”  T,  taken  to  be  0.2  (Gregg  1987).  Lien  et  al. 
(1998)  and  Lien  and  D’Asaro  (2006)  describe  a  method  to 
compute  £  from  the  Lagrangian  spectra  of  vertical 
acceleration  0a(ffl)  measured  by  a  Lagrangian  float.  Here, 
however,  the  high  noise  level  of  the  pressure  sensor 
prevents  this  method  from  being  used  except  in  a  very 
few  cases. 

Lor  very  weak  mixing  this  method  will  fail  as  the  inertial 
subrange  disappears.  The  numerical  study  of  D’Asaro  et  al. 
(2004)  suggests  that  this  limit  is  reached  at  values  of  £  below 
10“®  W  kg“',  equivalent  here  to  values  of  below  about 
10-6  m2  s-i. 


A  consistency  test — Using  from  Eq.  18  in  Eq.  12 

yields  the  simple  expression 

=  (23) 

If  mixing  is  entirely  diapycnal  Eq.  23  should  be  true. 
Deviations  indicate  either  a  significant  isopycnal  flux 
divergence  or  problems  with  either  the  budget  or  micro¬ 
structure  methods.  Using  Eq.  23  is  a  better  test  for 
diapycnal  balance  than  directly  comparing  with  Vn 
because  it  does  not  require  dividing  by  a  quantity  that 
can  be  zero,  and  because  it  eliminates  N  from  the 
calculation. 


%(ffl)  =  PsX  (19) 

It  is  white  and  proportional  to  the  dissipation  rate  for  the 
scalar.  The  value  of  the  Kolmogorov  constant  /?.,  =  0.6  has 
an  uncertainty  of  less  than  a  factor  of  two.  In  principle, 
then,  can  be  estimated  directly  from  the  level  of  the 
spectrum  of  the  Lagrangian  rate  of  change  of  potential 
density  0^  within  the  inertial  subrange. 

Practically,  the  finite  size  of  real  floats  leads  to 
deviations  of  the  spectral  forms  from  Eq.  19.  Observed 


Data  and  analysis 

The  Lagrangian  float — The  analysis  will  use  data  from  a 
single  Lagrangian  float,  MLEII  (mixed-layer  float,  second 
generation)  number  8  since  this  was  the  only  float 
measuring  chlorophyll.  The  float  profiled  every  12  h  from 
a  depth  of  40-60  m  to  the  surface.  Between  profiles,  i.e., 
during  “drifts,”  it  actively  adjusted  its  buoyancy  to  remain 
near  the  25.5  kg  m“3  isopycnal.  The  float  position  was 
determined  at  the  top  of  some  profiles  using  system 
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ARGOS  because  the  global  positioning  system  (GPS)  did 
not  function.  The  float  measured  pressure,  temperature, 
salinity  (Seabird  41CT),  downwelling  irradiance  at  490  nm 
(Biospherical  Instruments  QCP200),  and  chlorophyll  fluo¬ 
rescence  (Wetlabs  ECO)  every  100  s  during  drifts  and  with 
roughly  1-m  resolution  during  profiles.  Details  on  these 
floats  can  be  found  in  D’Asaro  (2003).  Details  on  their 
measurements  of  the  subinertial  and  internal  wave  envi¬ 
ronments  on  the  Oregon  shelf  can  be  found  in  D’Asaro 
(2004)  and  D’Asaro  et  al.  (2007). 

The  float  was  deployed  on  23  July  2001  (day  205.62)  at 
about  45.5°N  at  the  100-m  isobath  (Fig.  1).  The  float 
traveled  south  in  the  prevailing  current,  moving  offshore 
and  stalling  near  the  shelf  break  on  days  210-215,  then 
moving  more  rapidly  across  Heceta  Bank  by  day  220. 
During  this  time  the  target  isopycnal  remained  at  20^0  m. 
On  about  day  230,  the  float  moved  offshore  near  Cape 
Blanco  and  descended,  with  the  isopycnal,  to  a  depth  of 
about  50  m.  The  mission  continued  through  day  275,  but 
this  latter,  offshore  part  of  the  trajectory  is  not  considered 
here. 

During  the  “drift”  periods  between  profiles,  these  floats 
are  Lagrangian  at  high  frequencies  and  isopycnal  at  low 
frequencies  (D’Asaro  2003).  At  low  frequencies  they 
dynamically  adjust  their  volume  in  response  to  the 
measured  T,  S,  and  P  so  that  the  float  measures  a  constant 
potential  density  and  therefore  remains  on  an  isopycnal 
surface.  At  high  frequencies  they  are  Lagrangian,  i.e.,  they 
follow  the  three-dimensional  motion  of  the  water.  The 
transitional  frequency  between  Lagrangian  and  isopycnal 
behavior  is  about  N/30;  the  additional  drag  of  the  large 
horizontal  cloth  drogue  decreases  this  frequency  to  well 
below  N. 

Bio-optical  models — Scope:  This  section  describes  the 
calibration  of  the  irradiance  and  chlorophyll  fluorescence 
sensors  and  the  models  used  with  these  data  to  estimate  the 
rate  at  which  solar  radiation  heats  the  water. 

Radiometer  calibration:  The  Biospherical  QCP200  mea¬ 
sured  downwelling  irradiance  at  490  nm.  This  was  checked 
by  pointing  the  sensor  upward  at  a  clear  blue  sky  on  1939Z 
11  July  2001.  Using  the  manufacturer’s  calibration,  the 
sensor  measured  £490  =  1.64  W  nm“h  At  the  same 
time  a  shortwave  radiometer  about  1  km  away,  corrected 
for  aging  (D’Asaro  2007),  measured  a  total  shortwave 
radiation  =  911  W  m“'.  The  online  version  of  the 
coupled  ocean  and  atmosphere  radiative  transfer  (http:// 
snowdog.larc.nasa.gov/jin/rtset.html)  (Jin  et  al.  2006)  was 
used  with  the  following  parameters:  visibility  140  km  (from 
Puget  Sound  Clean  Air  Agency),  integrated  ozone  amount 
=  325  Dobson  units  (DU)  (from  Environment  Canada, 
ftp://es-ee.tor.ec.gc.ca/pub/ozone_maps/nh/),  a  “mid-lati¬ 
tude  summer”  atmospheric  model,  a  “MODTRAN  Mar¬ 
itime”  or  “OP AC  Continental  Clean”  mixed  layer  model, 
no  clouds,  and  a  total  precipitable  water  =  3.15  g  cm“2 
tuned  to  yield  the  observed  This  yielded  £490  = 
1.64  W  m-^,  about  3%  smaller  than  that  measured.  This  is 
within  the  expected  errors  in  radiometer  calibration. 


Fig.  1.  Track  of  float  8  (dashed  line)  overlaid  on  topography 
of  Oregon  shelf  Isobaths  at  50,  100,  150,  and  200  m  are  shown. 
Time  is  labeled  along  the  track;  every  10th  day  of  the  year  2001 
is  shown. 

Diffuse  attenuation  profile:  Vertical  profiles  of  irradi¬ 
ance  £490  were  made  twice  daily,  near  5Z  and  17Z.  The  5Z 
profiles  occurred  nearly  at  sunset  for  the  first  part  of  the 
mission.  Later,  they  occurred  after  dark  and  were  not 
useful.  The  17Z  profiles  occurred  during  mid-morning  with 
a  solar  elevation  near  30°.  As  can  be  seen  in  Fig.  2,  a 
complete  profile  often  required  using  both  the  up  and  down 
casts.  Typically,  however,  the  irradiance  from  the  up  and 
down  profiles  did  not  match,  presumably  because  the 
surface  irradiance  changed  during  the  approximately  1-h 
interval  between  them.  Because  the  shapes  were  usually 
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Fig.  2.  Two  typical  irradiance  profiles.  In  each,  the  upward- 
going  segment  (red)  is  offset  from  downward-going  segments 
(blue,  dotted  line)  occurring  about  1  h  later.  These  are  merged 
into  a  single  profile  by  multiplying  the  downward  segment  by  an 
appropriate  factor  (blue  solid  line).  This  works  well  in  case  B,  less 
well  for  case  A.  The  black  dashed  line  shows  the  fit  of  Eq.  24  to 
the  combined  profile  and,  for  case  A  only,  an  extension  to  very 
low  light  level  using  Eq.  25.  The  profile  for  case  B  is  shifted  by  a 
factor  of  10  so  that  it  does  not  overlap  the  profile  from  case  A. 
The  irradiance  level  equivalent  to  one  count  is  indicated.  The 
length  of  the  vertical  axis  on  this  plot  is  400  kPa,  which  is 
approximately  equal  to  40  m. 


very  similar,  the  two  profiles  were  combined  by  shifting  the 
down  profile  to  match  the  up  profile  using  a  factor 
determined  for  each  profile  by  eye.  Only  profiles  with  at 
least  10  data  points,  a  maximum  depth  of  at  least  5  dbar, 
and  a  maximum  irradiance  of  at  least  50  counts  were  used. 
Data  points  with  less  than  three  counts  were  excluded  from 
the  analysis.  The  combined  profile  was  fit  with  a  log- 
quadratic  function 

£490  =  (24) 

and  with  a  linear  form 

£490  =  £oe-"^^  (25) 

A  complete  irradiance  profile  was  constructed  using  Eq.  24 
at  all  depths  where  this  form  decreased  with  depth, 
occasionally  extending  the  profile  below  this  depth  using 
the  constant  value  of  £490  from  Eq.  25. 

Two  sample  profiles  and  fits  are  shown  in  Fig.  2.  The 
combined  profiles  are  accurately  fit  by  Eq.  24  for  case  B. 
For  case  A,  however,  £490  is  sufficiently  large  that  the  light 
level  falls  below  the  sensitivity  of  the  instrument  within  the 
profile.  This,  and  the  variability  of  the  surface  illumination 
with  time,  results  in  a  cusp  in  the  irradiance  curve  at  depth. 
This  could  probably  be  eliminated  by  more  careful  analysis, 
but  the  light  levels  are  so  low  that  these  details  have  little 
effect  on  water  column  heating,  the  goal  of  this  analysis. 


Diffuse  attenuation  coefficient  (m ') 


Fig.  3.  Scatter  plot  of  £490  and  calibrated  chlorophyll  high 
light  level  (dots).  Blue  line  is  the  Morel  and  Maritorena  (2001) 
relationship.  Red  line  shows  the  values  of  £490  averaged  in 
chlorophyll  bins. 

Fluorometer  calibration:  The  Wetlabs  ECO  fluorescence 
reading  £  was  converted  to  chlorophyll  using  the  Morel 
and  Maritorena  (2001)  empirical  relationship  for  type  1 
waters  between  chlorophyll  concentration  Chi  [mg  m“3] 
and  the  diffuse  attenuation  coefficient  at  490  nm  £490  [m“ '] 

£,„„,  =  0.01660  -t-  0.07242C/r/°“‘'^^  (26) 

This  method  is  insensitive  to  changes  in  the  irradiance 
sensor  calibration  due  to  biofouling  because  £490  depends 
only  on  the  ratios  of  irradiance  with  depth,  not  the  absolute 
irradiance. 

The  fluorescence  data  show  large  increases  beginning 
about  day  235.  The  increases  were  rarely  found  during 
down  profiles,  but  were  common  during  up  profiles, 
suggesting  that  they  were  due  to  a  piece  of  plant  attached 
to  the  float  that  was  swept  upward,  out  of  view  of  the 
fluorometer,  during  down  profiles.  The  increases  were 
removed  from  the  profile  data  by  deleting  all  data  from  up 
profiles,  which  was  more  than  a  factor  of  two  larger  than 
the  data  from  down  profiles  interpolated  to  its  location. 
The  analysis  shown  here  ends  on  day  240,  so  these  artifacts 
have  little  effect. 

Figure  3  shows  the  calibrated  Chi,  computed  using  Chi 
=  a  +  bF,  plotted  against  Chl^^.  The  constants  a  and  b 
were  adjusted  to  provide  a  good  fit.  Only  data  for  which 
£490  >  8  X  10“3  Wm“2  nm“'  (symbols)  are  plotted  in 
Fig.  3  and  used  as  data  because  smaller  values  do  not 
correlate  well.  The  average  accuracy  of  the  calibration, 
given  by  the  standard  deviation  of  <£>/£,„,„,  where  <£>  is 
£490  averaged  in  Chi  bins  (white  circles)  and  £,„,„  is  the 
Morel  and  Maritorena  (2001)  value  of  £490  for  each  bin,  is 
11%. 

Radiation  fields  and  heating  rate:  The  model  irradiance 
profiles  (e.g.,  Fig.  2)  were  interpolated  to  compute  the  ratio 
-E^49o(0)/£49o(-Z/)  of  irradiance  at  the  surface  and  at  the  float 
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Fig.  4.  Isopycnal  analysis  for  drift  segment  12,  a  segment  with  an  especially  large  diapyenal  correction.  Profiles  of  (a)  potential 
density  a,  (c)  chlorophyll,  and  (f)  spice  n  from  before  (black)  and  after  (green)  this  drift  segment.  Red  circles  show  data  during  drift. 
Time  series  during  drift  of  (b)  pressure  (red  solid  line,  400  kPa  =  40  dbar)  and  a  (black  solid  line),  (d)  chlorophyll  (green),  and  (g)  fl 
(green).  Scatter  plots  of  (e)  Chi  and  (h)  If  against  a  during  the  drift.  A  least-squares  quadratic  fit  (blue  line)  is  computed  after  recursively 
rejecting  data  deviating  by  more  than  3  SD  from  fit.  Rejected  data  are  plotted  in  black.  These  quadratic  functions  are  used  to  predict 
Chl(a  [t])  and  IT  (a  [?])  (magenta)  in  panels  d  and  g,  respectively,  and  to  compute  their  values  at  the  target  isopycnal  (blue)  C/i/q  and  Ho  in 
panels  d  and  g,  respectively.  The  rate  of  change  of  CA/q  and  ITo  is  computed  from  a  least-squares  line  (blue  dashed  lines  in  panels  d  and  g, 
respectively).  Target  isopycnal  co  is  shown  in  panels  b,  e,  and  h  (black  dashed  line). 


depth.  This,  times  the  measured  £'49o(Zy),  yields  a  nearly 
continuous  estimate  of  £'490(0)  except  at  times  when 
£49o(Zy)  is  so  small  that  reliable  estimates  cannot  be  made. 
Because  the  heating  at  these  times  is  also  very  small,  the 
large  relative  errors  are  not  important. 

The  full  visible  spectrum  (200-700  nm)  of  downward 
irradiance  at  the  surface,  £^(0),  was  computed  from  £490(0) 
using  the  spectral  shape  shown  in  fig.  74  of  Jerlov  (1976). 
The  shortwave  irradiance  above  the  surface  is  estimated 
from  the  integral  of  £^(0),  with  the  addition  of  a  60% 
infrared  component,  i.e.,  Es„  =  J  £a(0)J/1/0.6.  This  time 
series  was  compared  with  the  data  from  shortwave 
radiometers  on  a  mooring  (Boyd  et  al.  2002)  on  the  central 
Oregon  shelf  (44°59.757'N,  124°6.998'W).  Two  radiome¬ 
ters  on  the  mooring  differed  by  5%;  one  was  corrected  to 
the  other.  The  float  time  series  tracks  well  both  the  daily 
cycle  of  irradiance  and  the  variability  between  days  as  mea¬ 
sured  on  the  mooring.  Obvious  errors  occur  when  the  float 
is  very  dark,  with  some  differences  undoubtedly  related  to 
the  geographic  separation  of  the  float  and  mooring. 

During  days  222-240  £49o(Zy)  was  above  the  noise  level; 
the  average  values  of  £5,,,  from  the  mooring  and  the  float 
differed  by  only  5%  despite  being  separated  by  up  to 


250  km.  This  suggests  that  the  float-based  estimate  of  is 
accurate  to  5-10%. 

The  surface  spectrum  £>.(0)  was  propagated  downward 
using  the  Morel  and  Antoine  (1994)  empirical  expressions 
for  Kx  as  a  function  of  Chi.  The  total  downgoing  irradiance 
Ed{t,z)  was  computed  from  the  integral  of  the  spectrum  and 
the  local  heating  rate  Q{t,z)  from  its  change  with  depth. 
The  heating  at  490  nm  249o(6z)  was  also  computed.  Their 
ratio  £490  =  Q  /  Q490  was  computed  as  a  function  of  depth 
and  time  with  a  resolution  of  about  1  d.  A  time  series 
Q49o(Z/)  =  £49o(Zy)£490ram(C^/  [Zy])  was  measured  at  the 
float  at  much  higher  time  resolution.  The  heating  rate  at 
the  float  at  high  resolution  was  thus  estimated  as  Q{t)  = 

Q49o(Z/[f])R49o(6Zy). 

Mixing  budgets 

Methods — The  analysis  here  is  focused  on  evaluating  the 
terms  in  Eqs.  12,  23,  and  17  along  the  oq  =  25.5  kg  m“3 
isopycnal.  Analysis  was  conducted  separately  for  each 
approximately  9.8-h-long  (35.4  ks)  “drift.”  An  example  is 
shown  in  Fig.  4.  The  float  moved  vertically  about  20  m  and 
diapycnally  about  0.15kgm“3  during  this  segment 
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(Fig.  4a, b).  The  measured  chlorophyll  and  spice  time  series 
(Fig.  4d,g,  green)  show  variations  that  are  in  phase  with  the 
a  fluctuations  as  shown  by  the  tight  lines  in  scatter  plots 
(Fig.  4e,h)  of  Chi  and  Ft  against  a.  Least-squares  quadratic 
fits  were  made  to  these  relationships  (blue  lines).  For  spice, 
let  n^(o')  be  the  least-squares  quadratic  fit  from  Fig.  4h. 
This  predicts  a  variation  nz)(0  =  n,/o'[?])  due  to  the 
diapycnal  motion  of  the  float  alone.  This  is  plotted  as  the 
magenta  line  in  Fig.  4g,  with  a  corresponding  Chi  plot  in 
Fig.  4d.  As  expected  n£)(t)  and  Chlr>{t)  account  for  much 
of  the  measured  variability.  This  is  removed  by  computing 
no(0  =  n(t)  -  14/3(0  +  rig(o'o),  the  blue  line  in  Fig.  4g,  with 
a  corresponding  C/i/o(0  in  Fig-  4d.  These  have  much  less 
variability.  They  will  be  interpreted  as  a  time  series  of  IT 
and  Chi  at  the  uo  isopycnal. 

The  average  rate  of  change  on  the  isopycnal,  i.e.,  n,  and 
C,  in  Eqs.  12  and  17,  was  computed  for  each  segment  from 
the  slope  of  a  least-squares  fit  line  to  FlofO  and  Chlo{t), 
respectively.  The  error  in  this  estimate  was  computed  by 
assuming  that  the  data  are  the  sum  of  this  line  and  random 
correlated  noise  with  a  correlation  function  similar  to  the 
average  correlation  function  estimated  from  the  data. 
Random  realizations  of  this  model  were  computed  and 
confidence  limits  computed  from  the  distribution  of  the 
slope  of  the  lines  fit  to  these  realizations.  Notice  that 
AilT  =  Yltdt,  the  integrated  changes  on  the  isopycnal,  is 
not  the  same  as  A/rlT  =  —  n(fi),  the  total  change 

along  the  float  trajectory,  because  the  float  moves 
horizontally  during  the  profiling  and  makes  surface 
communications  between  each  drift.  Thus  each  drift  can 
be  considered  a  separate  9.8-h-long  Lagrangian  experiment 
starting  at  a  different  horizontal  location  on  the  isopycnal. 
A/,n  measures  the  cumulative  changes  from  these  experi¬ 
ments,  ignoring  the  changes  due  to  the  location  jumps 
between  experiments.  D’Asaro  (2004)  discusses  the  magni¬ 
tude  of  these  jumps  in  more  detail. 

The  quantities  11^^  and  are  computed  similarly.  For 
each  segment,  all  data,  both  drifts  and  profiles,  extending 
from  the  center  of  the  preceding  segment  to  the  center  of 
the  following  segment  and  within  6a  of  (Tq  are  fit  with  a 
fourth-order  polynomial  in  a  and  the  derivatives  evaluated 
at  (To-  As  will  be  discussed  below,  three  values  of  6a  from 
0.2  to  0.3  kg  m“3  were  used;  smaller  values  lead  to  a  large 
uncertainty,  and  include  a  domain  that  is  too  far  away  from 
the  target  isopycnal.  The  error  in  is  large  and  difficult 
to  estimate  as  there  is  no  good  understanding  of  the 
statistics  of  meter-scale  n  variations.  An  error  was 
estimated  by  making  separate  polynomial  fits  for  the  first 
and  second  halves  of  each  segment.  The  standard  error  was 
taken  as  the  absolute  value  of  the  difference  between  these 
two  estimates. 

The  stratification  N  for  each  float  drift  is  taken  from 
D’Asaro  et  al.  (2007)  (their  fig.  4a).  The  float  has  only  a 
single  conductivity-temperature-depth  (CTD)  sensor  so 
continuous  measurements  of  stratification  are  not  avail¬ 
able.  Values  of  N  taken  from  the  twice-daily  profiles,  e.g., 
Fig.  4a,  are  highly  variable  and  do  not  yield  accurate 
results.  D’Asaro  et  al.  (2007)  instead  use  the  high-frequency 
cutoff  of  the  internal  wave  spectrum  to  compute  a  robust 
average  N  for  each  drift  segment.  The  typical  difference 


between  this  estimate  and  the  CTD-based  estimate  for  the 
one  float  for  which  this  is  available  is  30%. 

The  spice  budget — Figure  5  shows  the  diapycnal  mixing 
budget  for  If.  The  float  drifted  at  10-30-m  depth  (Fig.  5a) 
until  day  234,  when  it  dove  to  40-60-m  depth.  Within  each 
drift,  its  depth  fluctuated  over  a  10-20-m  range  because  of 
internal  waves.  However,  relative  to  the  target  isopyncal  (Tq 
=  25.5  kg  m“3  (Fig.  5b),  the  fluctuations  are  much  less, 
only  a  few  meters.  The  profiles  (red  curves  in  Fig.  5a)  map 
out  the  variations  in  H  surrounding  (Tq  as  shown  in  Fig.  5b. 
The  target  isopycnal  is  near  a  minimum  in  H,  a  cold,  fresh 
layer  most  likely  a  remnant  of  the  deep  winter  mixed  layer. 
The  layer  persisted  until  the  float  passed  into  a  different 
water  mass  on  day  234.  During  this  time  H^  is  positive, 
consistent  with  diffusive  warming  of  the  layer  (Fig.  5d).  A 
small  correction  for  solar  heating  (pink  lines)  makes  little 
change.  The  strong  curvature  of  H  near  its  minimum  results 
in  a  positive  value  of  (Fig.  5e).  The  dissipation  rate  of 
potential  density  variance  %  and  the  diffusivity 
computed  from  it  (Fig.  5f)  vary  over  two  orders  of 
magnitude.  Equation  23  predicts  similarly  large  variations 
in  the  rate  H,  at  which  diapycnal  diffusion  erases  the  cold 
fresh  intrusion  (Fig.  5g).  The  measured  H,  mostly  tracks 
these  variations  well  within  the  uncertainties.  Thus  a 
diapycnal  mixing  balance  applies  most  of  the  time.  A 
detailed  discussion  of  the  comparison  is  deferred  to  the  next 
section. 

This  calculation  depends  on  the  choices  of  time  and 
space  averaging.  Given  the  float  sampling  scheme,  the 
natural  timescale  is  the  drift  duration.  However,  averaging 
over  about  four  such  intervals,  about  1.8  d,  is  necessary  to 
reduce  the  noise.  For  a  diffusivity  of  10“4  (10“^)  m^  s“', 
this  corresponds  to  a  diapycnal  scale  of  4  m  (3  m)  and, 
for  N  =  0.015  s“',  a  density  difference  of  Acr  = 
0.1(0.07)  kg  m“3.  However,  values  of  H^^  computed  over 
density  ranges  this  small  are  quite  noisy;  the  larger  values  of 
Aa  used  in  Fig.  5  produce  more  stable  results,  but  may  also 
introduce  a  bias  in  the  value  of  H^^. 

The  relatively  good  match  between  the  two  estimates  of 
H^  in  Fig.  5g  was  achieved  only  by  increasing  the  value  of 

to  1.0  from  the  D’Asaro  and  Lien  (2007)  value  of  0.6. 
This  increased  the  estimated  value  of by  a  factor  of  1.6. 
The  change  is  well  within  the  current  uncertainty  in  and 
could  represent  a  more  accurate  estimate  of  its  value. 
However,  it  could  also  be  due  to  a  low  bias  in  the  estimate 
of  H^^. 

The  effect  of  using  a  nonlinear  and  thus  nonconservative 
H  in  these  calculations  was  assessed  by  defining  a  new  spice 
function  for  each  drift  segment  of  the  analysis,  linearized 
around  the  mean  values  So  and  0o  for  that  segment, 

n{s,e) = n(So,6io) + ns(So,eo)  [s  -  So]  +  +  n0(So,eo)  [o  -  do]  (27) 

where  H^  and  Hq  are  the  partial  derivatives  of  H.  Figure  5g 
shows  the  results  of  using  ri(S,0)  instead  of  n(S,0)  using  a 
dimmed  version  of  the  same  lines.  These  two  different 
computations  are  barely  distinguishable,  showing  that  the 
use  of  a  conservative  H  has  little  effect. 


j  j  Diapycnal  Pressre 

(10‘"kgm'^s‘')  D^^andZ  m^kg'*  n  (10‘^  kg  m'^  s"’)  n'*  (10'^  s"'’)  displacement  (m)  (kPa) 
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An  intense  mixing  event — A  striking  feature  in  the  data  is 
the  intense  mixing  event  on  day  222,  examined  in  detail  and 

da 

at  maximum  resolution  in  Fig.  6.  Spectra  of  —  (Fig.  6a) 

show  a  spectral  level  during  the  drift  segment  centered  on 
day  222.4  about  25  times  larger  than  the  average  value  in 
the  four  nearby  segments.  Assuming  /i.,  =  0.6  and  F  =  0.2, 
this  spectrum  implies  e  =  9  X  lO-’^m^s"^,  a  diapycnal 
diffusivity  =  9  X  10“4  m^  s“i  averaged  over  the  T  = 
10.8-h-long  segment,  and  a  mixing  thickness  5  of 

about  6  m.  A  partial  homogenization  of  the  water  column 
over  a  depth  of  5-10  m  is  seen  by  the  reduced  N  (Fig.  5c) 
and  the  reduced  minimum  in  Ft  (Figs.  5b,  6e)  extending  for 
about  1  d  after  the  event. 

This  mixing  is  intense  enough  to  allow  a  second  estimate 
of  the  mixing  rate  to  be  made  using  the  acceleration 
spectrum  (Fig.  6b).  The  spectra  for  day  222.4  clearly  stand 
out  above  the  adjoining  spectra  for  oAN  >  2.  For  each  drift, 
the  acceleration  spectrum  is  predicted  from  e  assuming  a 
Froude  number  coq/N  =  0.5  (Lien  and  D’Asaro  2006)  and 
plotted  in  Fig.  6b  as  a  dashed  line.  The  predicted  spectrum 
matches  the  measured  spectrum  well  on  day  222.4, 
implying  that  the  two  independent  diffusivity  estimates 
Vs  (Eq.  18)  and  (Eq.  22)  are  nearly  the  same.  At  the 
other  times,  the  measured  spectrum  is  dominated  by 
pressure  sensor  noise  (black  dashed  line)  and  internal 
waves,  with  the  predicted  turbulent  signal  falling  far  below 
these.  Only  the  V^  estimate  is  useful. 

Figure  6c,d  shows  the  pressure,  a,  and  Tt  variations  for  the 
three  drift  segments  centered  on  the  event.  Before  and  after 
the  event,  the  float  remains  witliin  0.03  kg  m“3  of  the  target 
isopycnal  and  FI  varies  by  less  than  0.01  kg  m“3  from  its 
average  value  in  each  drift.  At  about  day  222.31  the  measured 
density  rapidly  decreases  by  about  0.1  kg  m“3^  corresponding 
to  a  diapycnal  displacement  of  about  5  m  upward  (Fig.  5b) 
and  then  returns  to  the  target  isopycnal  over  the  next  2.5  h. 
The  value  of  n  similarly  increases,  but  does  not  return  to  its 
initial  value,  instead  settling  at  a  value  about  0.04  kg  m“3 
higher.  This  results  in  an  obviously  positive  value  of  IT, 
shown,  for  this  drift  segment,  by  the  dashed  line. 

Figure  6e  shows  n((T)  and  the  polynomial  fit  over 
0.3  kg  m“3  used  to  estimate  Ft,,,,.  Figure  6f  shows  the 
resulting  budget  of  Ft,,  as  in  Fig.  5g.  The  measured  value  of 
n,  (black)  agrees  with  that  estimated  (colors)  from  Eq.  23 
within  the  errors. 


These  data  are  interpreted  as  follows:  The  water  parcel 
tracked  by  the  float  is  entrained  into  the  bottom  of  a  5-10- 
m  thick  mixing  patch  at  about  day  222.31.  The  simulation 
of  D’Asaro  et  al.  (2004),  their  fig.  3c,  shows  the  Lagrangian 
view  of  a  similar  event.  Because  the  mixing  patch  lies  above 
the  float’s  initial  isopycnal  ao,  the  mixed  water  is  both 
lighter  and  has  a  higher  Tt  than  that  initially  at  ao.  Because 
the  float  is  programmed  to  be  isopycnal,  it  sinks  back  to  ao 
over  the  few  hours  after  entrainment,  but  remains  within 
the  patch  of  mixed,  higher-lT  water.  The  presence  of  an 
organized  packet  of  internal  waves  starting  at  day  222.4 
suggests  a  larger  structure  associated  with  the  mixing  event. 
Otherwise  no  explanation  of  its  cause  is  presented  here. 

The  chlorophyll  budget — Figure  7  shows  a  similar  anal¬ 
ysis  for  chlorophyll  using  the  diagnostic  equation 

Chl,=  '^Chh„  (28) 

Estimates  of  Chl^a  and  x  are  computed  as  described  above 
for  n.  The  observed  rate  of  change  on  an  isopycnal  (left  in 
equation,  black  in  Fig.  7e)  is  compared  with  that  expected 
from  diapycnal  diffusivity  alone  (right  in  equation,  colored 
in  Fig.  7e).  The  agreement  is  poor,  indicating  that  Eq.  28 
does  not  apply  and  that  factors  other  than  diapycnal 
mixing  are  important  in  changing  the  chlorophyll  concen¬ 
tration. 

Figure  7d  compares  the  chlorophyll  measured  by  the 
floats  (green  circles,  only  every  10th  plotted)  with  the 
cumulative  change  in  chloropyll  from  the  drift  segments 
only  (black).  For  example,  near  day  220  the  green  dots 
show  chlorophyll  levels  decreasing  past  5  mg  m“3^  corre¬ 
sponding  to  the  contoured  levels  in  Fig.  7a, b  at  the  float 
(yellow).  However,  the  trends  during  each  drift  are 
increasing  (red  lines)  so  that  the  cumulative  change  during 
drifts  (black  line,  ‘Lagrangian  Change’)  indicates  an 
increasing  chlorophyll  during  the  drifts  even  though  the 
measured  chlorophyll  is  overall  decreasing  with  time.  This 
difference  is  explained  by  the  large  decreases  that  occur 
between  drifts.  During  the  profdes  and  surfacing  the  float  is 
moved  horizontally  relative  to  the  target  isopycnal  so  when 
it  again  settles  onto  the  isopycnal,  it  is  at  a  different 
location  and  measures  a  different,  and  in  this  case 
decreased,  value  of  chlorophyll.  These  changes  do  not 
represent  changes  following  the  water,  rather  they  show  the 
presence  of  spatial  gradients.  It  is  the  Lagrangian  changes. 


Fig.  5.  Computation  of  the  diapycnal  budget  of  spice  IT.  (a)  Float  trajectory.  Drifts  are  black,  profiles  are  red.  (b)  Displacement  of 
float  from  target  isopycnal  cq  =  25.5  kg  m“3  during  drifts  (yellow),  map  of  spice  FI  around  cto  (color,  black  solid  line  contours),  position 
of  isopycnals  above  and  below  (Tq  (white  lines)  by  0.2,  0.3,  and  0.4  kg  m“3.  (c)  Stratification  at  float  computed  for  each  drift  (solid  line) 
and  averaged  over  about  1.8  d  (dashed  line),  (d)  Rate  of  change  of  spice  at  gq  averaged  as  in  panel  c  corrected  for  solar  heating  (black 
solid  line)  and  uncorrected  for  solar  heating  (pink  solid  line).  Vertical  bars  show  1  and  2  SD  of  error,  (e)  Second  derivative  of  spice  with 
respect  to  potential  density  computed  over  three  ranges  in  Act  as  indicated  by  the  line  colors.  Dashed  lines  are  computed  over  two 
drift  segments;  solid  lines  are  averaged  over  1.8  d.  Vertical  bars  show  1  and  2  SE.  (f )  Dissipation  rate  of  potential  density  variance  (black 
solid  line)  computed  using  Eq.  20.  Diapycnal  diffusivity  (blue  solid  line)  computed  from  x  using  Eq.  8.  Both  are  averaged  over  1.8  d. 
Vertical  lines  are  68%  and  95%  confidence  limits  of  the  fits,  (g)  Test  of  diapycnal  balance  using  Eq.  23.  FI,  from  panel  d  (black  solid  line) 
is  compared  with  xl’^l^aa  (colored  lines)  from  panels  e  and  f  adjusted  to  increase  the  value  of  [i,  to  1.0  from  the  nominal  0.6.  Vertical  bars 
are  1  SE.  All  plots  in  panel  g  were  also  computed  using  a  linearized  version  of  Ft  and  plotted  in  panel  g  as  light  lines  of  the  same  color 
slightly  offset  in  time.  This  is  most  apparent  on  the  red  lines  near  day  210. 
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Fig.  6.  Detailed  view  of  the  intense  mixing  event  on  day  222.  (a)  Spectra  of  potential  density  change  for  the  intense  mixing  drift 
(black)  and  surrounding  drifts  (colored).  Model  spectra  (dashed  line)  are  fit  for  co  <  2a)/,.  The  value  of  oj/,  is  indicated  by  the  heavy 
circles,  (b)  Spectra  of  vertical  acceleration  for  the  event  drift  (black)  and  surrounding  drifts  (colored).  Model  acceleration  spectra  (dashed 
lines)  with  e  computed  from  the  density  spectra,  (c)  Time  series  of  depth  and  a  for  the  event  drift  and  the  drifts  preceding  and  following  it. 
(d)  Time  series  of  spice  fl  for  these  three  drifts,  (e)  fl  as  a  function  of  a  for  these  three  drifts  (dots  and  circles)  and  quartic  fits  to  these 
(lines),  (f )  Spice  budget  for  the  3  d  centered  on  the  event  as  in  Fig.  5g. 
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a)  Chlorophyll  and  float  trajectory 
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Fig.  7.  The  diapycnal  budget  of  chlorophyll,  (a)  Color  map  and  contours  of  chlorophyll  Chi  from  profiles,  overlaid  by  the  float 
trajectory  (yellow)  shown  as  mean  and  standard  deviation  of  depth  for  each  drift  segment,  (b)  Color  map,  contours,  and  float  trajectory 
as  in  panel  a  but  relative  to  the  target  isopycnal  co  =  25.5  kg  m“3.  Positions  of  isopycnals  0.2,  0.3  and  0.4  kg  m“3  above  and  below  ffg 
(white  lines),  (c)  Second  derivative  of  chlorophyll  Chl„„  computed  over  three  ranges  of  Ac  as  indicated  by  the  line  colors  and  averaged 
over  1.8  d.  Vertical  bars  show  1  and  2  SE.  (d)  Chlorophyll  during  the  float  drifts  (green  dots,  only  every  10th  shown),  linear  fits  to  these  as 
in  Fig.  4  (red  lines).  Cumulative  Lagrangian  change  (black  line)  created  by  placing  the  red  lines  end  to  end.  (e)  Test  of  diapycnal  balance 
using  Eq.  28.  Rate  of  change  of  chlorophyll,  Chi,,  averaged  over  1.8  d  (black)  is  compared  with  (colored  lines)  similarly  averaged. 

The  value  of  has  been  increased  to  1.0  from  the  nominal  0.6.  Vertical  bars  are  1  SE. 
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the  cumulative  changes  during  the  float  drifts,  that 
represent  true  temporal  changes  following  the  water.  It  is 
these,  black  in  Fig.  7d,e,  that  do  not  agree  with  the 
predictions  of  Eq.  28  in  Fig.  7e. 


Interpretation 


Figure  5g  provides  a  quantitative  test  of  whether  spice, 
and  thus  implicitly  temperature  and  salinity,  are  governed 
by  the  diffusion  equation 


DYl 


+  G" 


(29) 


where  the  time  derivative  is  taken  in  a  water-following 
(Lagrangian)  frame  and  the  d  derivative  is  taken  perpen¬ 
dicular  to  the  isopycnal.  The  differences  between  Eqs.  29 
and  23  or  12  result  from  the  latter  two  being  set  in  the  more 
computationally  convenient  isopycnal  frame  rather  than 
Eangrangian  frame.  The  test  applies  for  measurements 
averaged  over  several  meters  and  several  days. 

If  Eq.  29  is  not  true,  there  are  several  possibilities.  Eirst, 
and  most  likely,  isopycnal  diffusion  is  important.  Other 
possibilities  include  different  diffusivities  for  temperature 
and  salinity  (McDougall  1984)  or  possibly  unknown  errors 
in  the  computation. 

Eigure  5g  shows  regions  where  diapycnal  mixing  dom¬ 
inates  and  regions  where  it  does  not.  The  clearest  diapycnal 
case  occurs  on  day  222  (the  first  vertical  yellow  bar)  and  is 
described  in  great  detail  above.  The  clearest  nondiapycnal 
region  is  before  day  210.  Here  H^^  is  clearly  positive  and  x 
is  large,  yet  H,  is  nearly  zero,  suggesting  that  the  tendency 
for  diapycnal  mixing  to  decrease  H  is  roughly  balanced  by 
other  processes,  presumably  isopycnal  mixing.  The  periods 
of  very  low  mixing  on  days  219  and  224—227  do  not  show  a 
diapycnal  balance.  However,  x  is  so  small  that  it  is  not  clear 
that  the  inertial  subrange  methods  are  accurate  here. 

The  most  interesting  period  is  associated  with  the  rapid 
increase  in  H  on  days  232-235  (the  second  vertical  yellow 
bar).  This  change  most  likely  represents  the  float  crossing  the 
front  between  subpolar  and  subtropical  waters.  A  front  is 
evident  in  satellite  sea  surface  temperature  images  near  the 
float’s  location  at  this  time.  However,  because  the  float’s 
location  is  poorly  known,  its  trajectory  relative  to  this  front 
cannot  be  determined.  The  float’s  isopycnal  rises  to  the 
surface  on  the  low  H  (cold)  side  of  the  front,  and  then 
plunges  to  about  50-m  depth  on  the  high  H  (warm)  side.  The 
float  itself  makes  this  transition  during  a  downward  profile, 
not  during  a  drift,  so  that  the  transition  is  not  Eagrangian 
and  does  not  necessarily  represent  the  path  by  which  water 
crosses  the  front.  As  the  float  approaches  the  front,  the  layer 
of  low  H  water  becomes  increasingly  thin  and  H^^  increases. 
The  mixing  balance  remains  nearly  diapycnal,  with  H,  also 
increasing  until  day  233.  As  the  float’s  isopycnal  surfaces 
and  plunges  switches  sign,  x  increases,  yet  H,  remains 
positive.  At  the  front,  therefore,  a  diapycnal  balance  no 
longer  holds,  despite  strong  diapycnal  mixing,  suggesting 
that  isopycnal  mixing  is  dominant  for  this  short  period  as 
also  reported  by  Alford  et  al.  (2005).  The  strong  mixing  also 
implies  a  strong  diapycnal  velocity,  estimated  at  a  few 


millimeters  per  second,  comparable  with  those  measured  by 
Barth  et  al.  (2004)  in  a  similar  environment. 

The  chlorophyll  budget  (Fig.  7e)  is  dominated  by  a 
bloom  extending  from  days  210  to  225.  The  float  skirted 
the  lower  edge  of  this  bloom.  As  discussed  above,  the 
Eagrangian  changes  within  each  drift  are  different  from  the 
changes  following  the  float.  The  Eagrangian  changes 
(Eig.  7d,  heavy  black  line)  show  an  increase  to  about 
7  mg  m“3  in  10  d  with  no  changes  thereafter.  The 
measured  chlorophyll  increases  irregularly  to  the  same 
value  and  then  decreases  to  almost  zero  in  about  5  d.  This 
latter  decrease  must  represent  a  spatial  gradient  in 
chlorophyll.  Detailed  SeaSoar  surveys  in  this  same  region 
(O’Malley  et  al.  2002)  just  after  the  float’s  transit  show  a 
similar  bloom,  with  its  outer  edge  corresponding  roughly  to 
the  isobath  of  the  float  (70-100  m).  The  gradient  at  the 
outer  edge  is  strong,  with  7  mg  m“3  change  occurring  in 
about  5  km.  Unfortunately,  the  failure  of  GPS  measure¬ 
ments  on  this  float  prevents  an  accurate  measurement  of 
the  surface  currents  and  thus  the  deviation  of  its  track  from 
the  water  motion  on  its  isopycnal.  However,  D’Asaro 
(2004)  estimates  drifts  for  float  10,  which  followed 
approximately  the  same  track,  as  roughly  1  km  d“i  during 
this  time.  It  is  thus  reasonable  that  the  float  traveled  just 
below  and  just  to  the  west  of  this  bloom,  in  a  region  of 
strong  isopycnal  and  diapycnal  gradients  in  chlorophyll. 

The  dominant  feature  of  the  Lagrangian  chlorophyll 
change  Chi,  (Eig.  7e)  is  the  increase  in  chlorophyll  as  the 
float  enters  the  bloom.  This  is  not  explained  by  diapycnal 
diffusion.  Could  it  be  due  to  phytoplankton  growth?  Light 
levels  at  the  float  are  low;  the  average  downward  irradiance 
from  days  210  to  220  is  about  2  W  m“2  at  the  float.  Using 
the  parameters  from  Spitz  et  al.  (2003)  for  light-limited 
growth,  this  would  correspond  to  an  e-folding  time  of  20  d 
compared  with  the  observed  growth  rate  {Chi,! Chi)  of 
about  4  d.  Growth  could  play  a  role,  particularly  if  the 
phytoplankton  at  this  depth  were  highly  light  adapted. 
However,  given  the  strong  correspondence  between  the 
observed  increase  at  the  float  and  the  presence  of  the  bloom 
above  the  float,  transport  is  a  more  likely  explanation. 
Sinking  of  a  few  meters  per  day  of  phytoplankton  could 
also  easily  explain  the  observed  increase  at  the  float.  The 
pattern  of  chlorophyll  on  days  213-218  suggests  this. 
Equally  likely  is  isopycnal  mixing  across  the  strong 
isopycnal  gradients.  Using  an  isopycnal  diffusivity  of 
4  m2  s“'  measured  in  this  area  (Dale  et  al.  2006)  and  a 
plausible  WChl  =  7  mg  m“3  /  (1,000  m)2  yields  Chi,  =  1.2 
X  10“^  mg  m“2  s“',  comparable  with  the  observations. 

Large  errors  in  the  second  derivatives  H^^  and  Chl„„  are 
major  weaknesses  in  this  analysis.  The  use  of  well- 
calibrated  pairs  of  CTD  sensors  on  the  floats,  as  commonly 
done  in  other  deployments,  would  help  significantly.  The 
other  major  weakness  is  the  lack  of  horizontal  measure¬ 
ments  to  define  isopycnal  processes.  Mapping  measure¬ 
ments  around  a  drifting  float,  similar  to  those  made  by 
Barth  et  al.  (2004),  is  the  most  obvious  approach. 

Data  from  an  isopycnal-following  float  were  used  to 
investigate  the  relative  importance  of  diapycnal  and 
isopycnal  mixing  in  changing  water  properties.  An  expres¬ 
sion  for  the  rate  of  change  of  a  scalar  measured  along  an 
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isopycnal  trajectory  due  to  diffusion  and  requiring  only 
diapycnal  derivatives  was  derived  following  McDougall 
(1984).  The  diapycnal  diffusivity  was  estimated  from  the 
spectrum  of  potential  density  at  frequencies  above  N 
(D’Asaro  and  Lien  2007).  Together,  these  form  a  test 
Eq.  23  for  a  purely  diapycnal  mixing  balance. 

This  test  was  applied  to  two  scalars,  spice  (Flament 
2002),  a  convenient  combination  of  temperature  and 
salinity,  and  chlorophyll.  Both  were  measured  for  35  d  by 
a  float  tracking  a  mid-depth  isopycnal  on  the  Oregon  shelf. 
Time  series  of  11  and  Chi  on  the  isopycnal  were  constructed 
by  correcting  for  the  small  diapycnal  motion  of  the  float. 
The  diapycnal  derivatives  were  estimated  from  the  float 
data  near  the  isopycnal.  A  bio-optical  model  was  used  to 
calibrate  chlorophyll  and  to  derive  the  local  rate  of  solar 
heating  from  measurements  of  downwelling  radiation  at 
490  nm  on  the  float.  Independent  tests  of  whether  the 
observed  rates  of  change  were  due  only  to  diapycnal 
diffusivity  were  generated  every  1.8  d.  The  following 
conclusions  about  mixing  apply  to  averages  of  these  data 
of  about  1.8  d  in  time  and  about  10  m  in  the  vertical. 

Variations  in  spice,  and  by  inference  temperature  and 
salinity,  were  primarily  due  to  diapycnal  mixing,  most 
spectacularly  in  a  single  large  mixing  event.  The  observed 
diapycnal  balance  could  only  be  made  to  work  on  average, 
however,  by  postulating  a  factor  of  1 .6  change  in  the  D’Asaro 
and  Lien  (2007)  estimate  of  the  appropriate  Kolomogorov 
constant.  Although  this  change  is  within  the  uncertainty  of 
their  estimate  and  perhaps  subject  to  additional  biases,  it 
emphasizes  the  continued  need  to  refine  such  estimates. 

At  times,  a  diapycnal  balance  did  not  apply  to  spice  and, 
by  inference,  isopycnal  mixing  must  have  been  important. 
Variations  in  chlorophyll  rarely  obeyed  a  diapycnal 
balance.  Instead,  some  unknown  combination  of  biological 
growth,  planktonic  sinking  across  the  strong  vertical 
gradients,  or  isopycnal  mixing  across  the  strong  isopycnal 
gradients  must  have  dominated.  The  observed  difference 
between  the  behavior  of  a  physical  and  biological  variable 
may  be  coincidence;  the  time  series  are  not  long  enough  to 
justify  a  more  general  conclusion. 

These  data  clearly  show  a  high  degree  of  variability  in 
diapycnal  mixing  in  the  thermocline  over  the  Oregon  shelf, 
even  when  averaged  over  several  days  (e.g.,  Fig.  5f).  This  is 
consistent  with  previous  measurements  (Gregg  et  al.  1993; 
Mourn  and  Nash  2000).  They  also  suggest  that  the  balance 
between  isopycnal  and  diapycnal  mixing  may  be  equally 
variable,  with  isopycnal  effects  perhaps  concentrated  at 
fronts,  and  with  the  relative  importance  of  isopycnal  and 
diapycnal  effects  for  scalar  mixing  highly  dependent  on  the 
distribution  of  the  scalar. 
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